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Abstract. Coupled dark matter-dark energy systems can suffer from non-adiabatic instabilities at early times and large scales. 
In these proceedings, we consider two parameterizations of the dark sector interaction. In the first one the energy-momentum 
transfer 4- vector is parallel to the dark matter 4- velocity and in the second one to the dark energy 4- velocity. In these cases, 
coupled models which suffer from non-adiabatic instabilities can be identified as a function of a generic coupling Q and of 
the dark energy equation of state w. In our analysis, we do not refer to any particular cosmic field. We confront then a viable 
class of models in which the interaction is directly proportional to the dark energy density and to the Hubble rate parameter to 
recent cosmological data. In that framework, we show that correlations between the dark coupling and several cosmological 
parameters allow for a larger neutrino mass than in uncoupled models. 
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INTRODUCTION 

Interactions between dark matter and dark energy are still allowed by observational data today. At the level of the 
background evolution equations, one can generally introduce a coupling between these two sectors as follows: 

Pdm + ^Pdm = Q, (1) 
p de + 3,^p de (l+w) = -Q. (2) 

Pdmipde) denotes the dark matter (dark energy) energy density, the dot indicates derivative with respect to conformal 
time dx = dt/a, Jti? = a/a and w = P de /p de is the dark-energy equation of state (P denotes the pressure). We 
work with the Friedman-Robertson- Walker (FRW) metric, assuming a flat universe and pressureless dark matter 
w dm = P dm j p dm = 0. 

Q encodes the dark coupling and drives the energy exchange between dark matter and dark energy. For e.g. Q < 
the energy flows from dark matter to dark energy. It also changes the dark matter and dark energy redshift dependence 
acting as an extra contribution to their effective equation of state. For e.g. Q < 0, dark matter redshifts faster, as a 
consequence, there is more dark matter in the past compared to uncoupled scenarios assuming that the dark matter 
density today is the same in the two models (see also the discussion in Ref. [1]). This general feature of coupled models 
is sketched in Fig.Q~](for a particular form of Q, see also Fig. [2J. 

Also notice that a universe in accelerated expansion today requires w < —1/3 even in the presence of a dark 
coupling. Indeed, the deceleration parameter satisfies 

4=-^2=\( l + 3w£l de) , (3) 

either with or without dark coupling. The curvature contribution has been neglected in this equation. 

In order to deduce the evolution of density and velocity perturbations in coupled models, we need an expression of 
the energy transfer in terms of the stress-energy tensor: 

V M ^ m)v =Gv and V,T$ e)v = -Q v . (4) 

The 4-vector Q v governs the energy-momentum transfer between the dark components and Tf~ dm \ v and are the 

energy-momentum tensors for the dark matter and dark energy, respectively. Equation (01) guaranties the conservation 
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FIGURE 1. This figure illustrates roughly the effect of a negative dark coupling term Q on dark matter and dark energy evolution, 
a is the scale factor, p,- with i = m,r,de corresponds to matter (dark matter plus baryon), radiation and dark energy densities 
respectively. Solid curve account for an uncoupled model with constant dark energy equation of state w > — 1 and dotted curve for 
coupled models with Q < with w > — 1. 



of the total energy-momentum tensor. In the following we consider two scenarios which reproduce at the background 
level Eqs. (Q~|i and (|2): 

1. Q v is parallel to the dark matter four velocity Uy ^: 

Qv = Qu { v dm) /a. (5) 

This choice of parameterization, which was first proposed in Ref. [2], avoids momentum transfer in the rest frame 
of dark matter (i.e v| dm j = 0), 

2. <2v is parallel to the dark energy four velocity Uy : 

Qv = Qu ( v e) /a. (6) 

In this case the dark matter velocity (Euler) equation is modified by the coupling and the weak equivalence 
principle is violated (see e.g. Ref [3] for a recent discussion). Notice that one example for this interaction is 

Q v = pPdin^v^ JUS, where (j) would be a coupled dark energy scalar field and °= V v 0. 

It was first pointed out in Ref. [2] that the dark coupling terms which appears in the non-adiabatic dark energy 
pressure perturbations are a source for early time instabilities at large scales for coupled models satisfying Eqs. 
and ©. Their analysis was however restricted to Q > proportional to the dark matter density, in which case the 
coupled model is particularly unstable for a constant dark energy equation of state. Non-adiabatic instabilities were 
subsequently analyzed by several authors, and it was proved that the stability depends on the type of dark coupling 
<2, on the dark energy equation of state w and on the Q v 4- velocity dependence, see Ref. fla LZL la l2D (for similar 
instabilities pointed out in coupled quintessence models, see Ref. OlQlllllD . 

Among these references, we proposed in Ref. [8], a criteria associated to the dubbed doom factor to identify the 
stability region of coupled models with a constant dark energy equation of state w satisfying Eqs. (0]l and (0. The 
doom factor is a function of the model parameters such as Q and w, but it is defined independently of the explicit form 
of the coupling Q. We review this result in the next section and extend it to models satisfying Eqs. (01 and (|6]). 

Based on this analysis, we study the compatibility of a successful class of models, in which Q is proportional to 
the dark energy density, with a fixed dataset including WMAP 5 year fUE! , HST Q, SN Hi, H(z) d] and 
LSS lfl7ll data. This work was done using the publicly available CAMB code lfl8tl and cosmomc package [19]. The 
latter were modified in order to include the interaction between the dark matter and dark energy components and the 
modified dark matter velocity equation in the case of coupled models satisfying Eq. (0 (see Eq. (fT8l). Present data 
will be shown to allow for a sizeable interaction strength and to imply weaker cosmological limits on neutrino masses 
with respect to non-interacting scenarios. 



ORIGIN OF NON-ADIABATIC INSTABILITIES 



Non-adiabatic instabilities arise at linear order in perturbations and appear to be driven by the dark coupling term 
present in the non-adiabatic dark energy pressure perturbation |2]. Using the definitions of gauge invariant density 
perturbation and entropy perturbation (see e.g. Ref. II20|P . one can work out a general expression which relates the 
dark energy pressure perturbation 8P £ i e in its rest frame to the one in any other frame. The former is characterized by 
^sde = [SPde/SPde] |rf> the propagation speed of pressure fluctuations in the rest frame of dark energy which has to be 
distinguished from c 2 ade = Pd e /Pde, the so called "adiabatic sound speed". In the synchronous or the Newtonian gauge, 
one obtains: 

Q 

§Pde = C 2 sde Sp de + (c 2 sde -C 2 ade )p de ^-, (7) 

where 8p c ( e denotes the dark energy density perturbation and 9 c / e = diV l r de \ is the divergence of the dark energy proper 
velocity, v l , de y Using equation (0, we see that the dark coupling resulting from the p de term directly affects the 8Pd e . 
In Eq. ([8]), we rewrite Eq. (Q in term of the dubbed doom factor d « Q which is a useful tool to spot the combined 
role of Q and w in driving the non-adiabatic instabilities. In the following we illustrate this feature in the framework 
of coupled models satisfying Eqs. <(5j and ©. 

Notice that non-adiabatic instabilities differ from the adiabatic ones. The latter appear at relatively small scales and 
late times. In the adiabati c reg ime, the effective sound speed of the fluid tends towards the adiabatic one which turns 
out to be negative I2TI I22, 23]. As a consequence, pressure no longer counteracts the effect of gravity and instabilities 
can breakout. 

In the following, we work with constant equation of state w in which case (c 2 d = w) and we assume that our 
universe is in accelerating expansion today, which implies that w < —1/3. Moreover, we restrict our analysis to the 
case c 2 de > and c 2 de = 1 will be assumed for numerical computation. 



GROWTH EQUATION AT LARGE SCALES AND THE DOOM FACTOR 

A cartoon equation of the growth equation governing the evolution of energy density linear perturbation for any species 
i,j is given by: 
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where 5,- = 5p,/p,- and the prime denotes a derivative with respect to the scale factor ' = d /da. The evolution of a 
perturbation depends on the relative weight of the three terms present in this equation and on their signs: 

1 . For positive A, the A and B terms taken by themselves would induce a rapid growth of the perturbation, which may 
be damped or antidamped (reinforced) depending on whether B is negative or positive, respectively. In particular, 
for A and B both positive, the solution may enter in an exponentially growing, unstable, regime. 

2. For negative A, in contrast, the A and B terms taken alone describe a harmonic oscillator, with oscillations damped 
(antidamped) if B is negative (positive). In the A,B < regime, the third term may play in fact the leading role. 

In the standard uncoupled scenario, the dark matter perturbations behave as in case Q] above (with A > and B < 0), 
while the dark energy ones provide an example of behavior as in case [2] 

For coupled models, we concentrate on the case in which the dark-coupling terms dominate over the usual one, in 
order to put forward the presence of non adiabatic instabilities (see Ref.(8|] for more details). For this purpose, let us 
rewrite Eq. (|7]i as: 

|^ = tf* + 3(«*-<S*)(l+w)(l+d)^, (8) 



where d refers to the doom factor which we have defined as: 



Q 

d= 3JPp de (l+w)' (9) 

The strong coupling regime can be characterized by |d| > 1, which also guarantees that the interaction among the two 
dark sectors drives the non-adiabatic contribution to the dark energy pressure perturbation. 

In Ref. [8], the strong coupling regime analysis was restricted to coupled models satisfying Eqs. and (|5). It is 
however rather straightforward to generalize this result to models satisfying Eq. (O (see also the appendix). At large 
scales-early times (J^/k ^> 1), the leading contributions in Q, or equivalently in d, to the second order differential 
equation for 8 de reads: 



SL - 3d(^+i)(% + 36%%fe^ + %^«[d])+... (10) 
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where we use the b notation introduced in Ref. |@1 , b — 1 stands for models with Q v °< u y ® and b = for models 



with Q v Uy le ^ ||6). The sign of the coefficient B e of 8 de in this expression is crucial for the analysis of instabilities. 



Assuming c 2 sde > 0, it reduces to the sign of the doom factor d. As previously argued, a positive d can trigger large scale 
instabilities. Similar second order differential equations for 8 de were obtained in Refs. laLZI] f° r particular expressions 
of the dark coupling Q and an analytical form of their solutions were derived in order to determine when 8 de blows 
up. In particular, the results of Ref. J3l confirm those of Ref. [2] for positive Q p dm and 1 + w > 0. In addition, 
our Eq. ( TTOb reproduces as well the dark coupling leading contributions of Ref. which studied Q °c p de coupled 
models, in particular, they carried out the first stability analysis in the framework of Q v oc p de u\? e ^ . 



VIABLE MODELS: Q °c p de 

Using the tools that we have developed in the previous section, we can now easily verify that for 

Q = $J?p de , (11) 

we have rather simple and viable models for specific combination of 1 + w and of the dimensionless constant coupling 
% . Indeed in these models the doom factor of Eq. © is given by: 

P 

(12) 



3(1 +w) 



When d < 0, that is, for | < and l+w>0(or^ >0 and 1 + w < 0), no instabilities are expected and we can safely 
fit the coupled models to cosmological data. Notice that this result is in agreement with those of Refs. 0, IB] which 
restricted though their stability analysis to the B, > case. 

For the sake of completeness we first give the expressions of the background dark fluids energy densities and first 
order perturbation evolution equations which were introduced in the CAMB code II 1811 . As previously mentioned, we 
work in the synchronous gauge. 



Background 

The solutions to Eqs. ([TJ and using Eq (fTTb are: 

Pdm = p d y-' + p^^(l~a-^)a-\ (13) 

Pde = p d °Ja-^ +w ^. (14) 

The dark energy density is thus always positive, all along the cosmic evolution and since its initial moment. To ensure 
that the same happens with the dark matter density, all values of w < are acceptable for £, < 0, while for positive £, 
it is required that 



In Fig. |2] we see that negative (positive) couplings lead to more (less) dark matter in the past than in the uncoupled 
case. In the following we focus on negative couplings in order to avoid non-adiabatic instabilities. 
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FIGURE 2. Scenario with Q °= p de . Relative energy densities of dark matter plus baryons £l d m+b (blue), radiation Sl rad (black) 
and dark energy Q. de (red), as a function of the scale factor a, for w=-0.9. Three values of the coupling are illustrated: £, = (solid 
curve), 0.25 (long dashed curve) and —0.25 (short dashed curve). 



Linear perturbation theory 

In the synchronous comoving gauge, metric scalar perturbations are described by the two usual fields [24] h(x, t) 
and T](x, t). Defining 8 = 8p/p for the fluid density perturbations, 9 = <9,v' for the divergence of the fluid proper 
velocity v' and using Eq. @, it results, at first order in perturbation theory: 
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(15) 
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(18) 



where b = 1 stands for models with Q v oc Uy* © and b — for models with Q v °= Uy ©. We have assumed that 
Jrff is the global expansion rate and that it does not contribute to Q perturbation. 

Notice that the continuity equation for dark matter <TT3T > always includes an extra term compared to uncoupled 
models which depends on the dark interaction Q of Eq. dTTb . Focusing on viable models, it can be showrQ that the dark 
matter growth equation at small scales is mainly modified in two ways: (i) the background evolution (in particular, p dm 



See 1 3] for similar coupled models. The details of the growth at small scales associated to the models studied here will be presented elsewhere. 
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FIGURE 3. Scenario with Q °c p^ e . The two upper panels correspond to models with Q v oc u v , the two lower panels to models 

with Q v oc u\f e \ Left ( right) panel: la and 2a marginalized contours in the ^-Q.^ m h^ ( £,-fv ) plane. The largest, green contours 
show the current constraints from WMAP ( 5 year data ), HST, SN and H (z) data. The smallest, red contours show the current 
constraints from WMAP (5 year data), HST, SN, H(z) and LSS data. 



and 3^ evolution) is different, (ii) the Hubble friction term and the source term get extra contributions from Q. For 
negative couplings, these two modifications lead to an enhancement of the dark matter growth compared to uncoupled 
models (see also Ref. [Q]] for similar models). The closer t, gets to -1, the larger is the growth. This particular feature 
can be constrained by cosmological data, in particular in the next section we will see that large scale structure data 
provide the strongest limits on the interaction Q. 

Also notice that the Euler equation for dark matter ( fTo*b is only modified in the Q v °= Uy^ 1 © case, leading a violation 
of the weak equivalence principle. Constraints resulting from the difference between dark matter and baryon velocities 
could provide additional restrictions on the allowed values of the dark coupling, to be added to the ones presented in 
the following section. For the sake of comparison with the already studied Q v °= Uy ^ © model we present below the 
likelihood plots for the two coupled models (© and ©. 

COSMOLOGICAL CONSTRAINTS FROM DATA FOR Q = %M>p de 

In Ref. dH], we explore the constraints on the dark energy-dark matter coupling £, using the publicly available package 
cosmomc QfJ. The latter is modified in order to include the coupling among the dark matter and dark energy 
components. More details on the cosmological model and on the priors adopted can be found in Ref. HI]. The datasets 



in the analysis are: 

1. WMAP 5-year data ifT^fll 

2. Prior on the Hubble parameter of 72 ± 8 km/s/Mpc from the Hubble key project (HST) Hill 

3. Super Novae la (SN la) data d 

4. H(z) data at < z < 1.8 from galaxy ages 

5. Large scale structure data (LSS data) from the Sloan Digital Sky Survey lfl7ll 

The data analysis is carried out into two runs, the first run includes the datasets from 1 to 4 while in the second run the 
fifth dataset is added. We restrict ourselves to w > — 1 and ^ < 0, a parameter region which ensures a negative doom 
factor, see Eq. (1121 1. and thus spans an instability-free region of scenarios to explore. 

Figure [3] (left panel) illustrates the 1 and 2a marginalized contours in the E,-Q. dm h 2 plane, where £t dm is today's 
ratio between dark matter energy density and critical energy density. The results from the two runs described above are 
shown. Notice that a huge degeneracy is present, being E, and £l dm h 2 positively correlated. The shape of the contours 
can be easily understood following our discussion in the previous sections. In a universe with a negative dark coupling 
<!;, there is an enhancement of the growth of structure relative to the non interacting case. The amount of intrinsic 
dark matter (which is directly proportional to Q. dm h 2 ) needed to reproduce the LSS data should decrease as the dark 
coupling becomes more and more negative. We also see that the addition of LSS data in the second run (red contours 
in Figs. [3]and[4]i gives the most stringent constraint on ^ . 

The right panel of Fig.[3]shows the correlation among the fraction of matter energy-density in the form of massive 
neutrinos f v and the dark coupling £, . The relation between the neutrino fraction used here / v and the neutrino mass 
for A^y degenerate neutrinos reads 

_ Q. v h 2 _ Lot v 1 _ N v m v 1 
/v ~~ a dm h 2 ~ 93.2eV ' £l dm h 2 ~ 93.2eV ' Q. dm h 2 ' ( 

Neutrinos can indeed play a relevant role in large scale structure formation and leave key signatures in several 
cosmological datasets. Degeneracies between dark energy sector parameters and £m v are rather well known, see 
Ref. ll25l l26l l27ll . Non-relativistic neutrinos in the recent Universe suppress the growth of matter density fluctuations 
and galaxy clustering. This effect can be compensated by the existence of a coupling between the dark sectors, due to 
the fact that in the coupled model negative couplings enhance the growth of matter density perturbations. 

Notice that the constraints on the Q v oc u v m) © and Q v °c Uy ^ (O models are rather equivalent, due to the fact that 
the background evolution history (see Eqs. (Tl3T > and (fl4b ) is the same, and the evolution of perturbations (fTBI l - dT8b is 
also very similar for these two models. 



CONCLUSIONS 

In these proceedings, we discuss the origin of non-adiabatic instabilities in coupled models satisfying 
^ '^dmv = Qu v mi " de ^ ja and V 1 pj dey — Qu v m ^ /a. We show that the sign of the doom factor 

d = ^^ n Q (^ \ (20) 
3JFp de (l +w) 

reveals the presence of instabilities. In particular, when d is positive and sizeable, d > 1, the dark-coupling dependent 
terms may dominate the evolution of the dark energy perturbation, and drive an unstable growth regime. Notice that d 
has been defined for a constant dark energy equation of state w (see e.g. Ref. |0,|2^] for w = w(a)). 

We have studied the constraints from current cosmological data on the dimensionless coupling ^ in a class of 
viable models in which Q — t,J^p de , independently of the dark interaction term 4-velocity dependence. The analyses 
presented here were carried out in the £, < and positive ( 1 + w) region of the parameter space, which offers the best 
agreement with data on large scale structure formation. From the results of our fits we find that both w and ^ are not 
very constrained from data: it can be noticed from Fig. [4] that substantial values for both parameters, near -0.5, are 
easily allowed. Furthermore, £, turns out to be positively correlated with Q. dm h 2 and a larger neutrino fraction f v than 
in uncoupled models is allowed for negative values of the coupling E, . 




FIGURE 4. Scenario with Q <x p de . Left (right) panel: la and 2a marginalized contours in the t,-w plane for Q v Uy" 1 ^ 

(Qv x My j. The largest, green contours show the current constraints from WMAP (5 year data), HST, SN and H(z) data. The 
smallest, red contours show the current constraints from WMAP (5 year data), HST, SN, H(z) and LSS data. 
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APPENDIX: GROWTH OF PERTURBATIONS IN STRONGLY COUPLED SCENARIOS 

The strong coupling regime can be characterized by 

° ' » |3(l+w)|, (21) 
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> 1-32 



sde\ ' 



(22) 



which ensure that the dark-coupling terms dominate the evolution of both 8 c i e and 9 c i e . With c 2 sde > 0, Eq. ( |2"TT i alone 
is enough to define the regime. 

For coupled models satisfying Eqs. ©, (|5]l and ©, the resulting growth equation for dark energy perturbations at 
large scales (Jf jk 3> 1) can be approximated by 
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where b = 1 for models with Q v « Uy' <(5j and b = for models with g v « u v e ' ©. 
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